\(\int x^8 (a+b x^4)^{3/4} \, dx\) [1032]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [A] (verified)
   Fricas [C] (verification not implemented)
   Sympy [C] (verification not implemented)
   Maxima [A] (verification not implemented)
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 15, antiderivative size = 125 \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=-\frac {5 a^2 x \left (a+b x^4\right )^{3/4}}{128 b^2}+\frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {5 a^3 \arctan \left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{256 b^{9/4}}+\frac {5 a^3 \text {arctanh}\left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{256 b^{9/4}} \]

[Out]

-5/128*a^2*x*(b*x^4+a)^(3/4)/b^2+1/32*a*x^5*(b*x^4+a)^(3/4)/b+1/12*x^9*(b*x^4+a)^(3/4)+5/256*a^3*arctan(b^(1/4
)*x/(b*x^4+a)^(1/4))/b^(9/4)+5/256*a^3*arctanh(b^(1/4)*x/(b*x^4+a)^(1/4))/b^(9/4)

Rubi [A] (verified)

Time = 0.03 (sec) , antiderivative size = 125, normalized size of antiderivative = 1.00, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.400, Rules used = {285, 327, 246, 218, 212, 209} \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=\frac {5 a^3 \arctan \left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{256 b^{9/4}}+\frac {5 a^3 \text {arctanh}\left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{256 b^{9/4}}-\frac {5 a^2 x \left (a+b x^4\right )^{3/4}}{128 b^2}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b} \]

[In]

Int[x^8*(a + b*x^4)^(3/4),x]

[Out]

(-5*a^2*x*(a + b*x^4)^(3/4))/(128*b^2) + (a*x^5*(a + b*x^4)^(3/4))/(32*b) + (x^9*(a + b*x^4)^(3/4))/12 + (5*a^
3*ArcTan[(b^(1/4)*x)/(a + b*x^4)^(1/4)])/(256*b^(9/4)) + (5*a^3*ArcTanh[(b^(1/4)*x)/(a + b*x^4)^(1/4)])/(256*b
^(9/4))

Rule 209

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[b, 2]))*ArcTan[Rt[b, 2]*(x/Rt[a, 2])], x] /;
 FreeQ[{a, b}, x] && PosQ[a/b] && (GtQ[a, 0] || GtQ[b, 0])

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 218

Int[((a_) + (b_.)*(x_)^4)^(-1), x_Symbol] :> With[{r = Numerator[Rt[-a/b, 2]], s = Denominator[Rt[-a/b, 2]]},
Dist[r/(2*a), Int[1/(r - s*x^2), x], x] + Dist[r/(2*a), Int[1/(r + s*x^2), x], x]] /; FreeQ[{a, b}, x] &&  !Gt
Q[a/b, 0]

Rule 246

Int[((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Dist[a^(p + 1/n), Subst[Int[1/(1 - b*x^n)^(p + 1/n + 1), x], x
, x/(a + b*x^n)^(1/n)], x] /; FreeQ[{a, b}, x] && IGtQ[n, 0] && LtQ[-1, p, 0] && NeQ[p, -2^(-1)] && IntegerQ[p
 + 1/n]

Rule 285

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[(c*x)^(m + 1)*((a + b*x^n)^p/(c*(m + n
*p + 1))), x] + Dist[a*n*(p/(m + n*p + 1)), Int[(c*x)^m*(a + b*x^n)^(p - 1), x], x] /; FreeQ[{a, b, c, m}, x]
&& IGtQ[n, 0] && GtQ[p, 0] && NeQ[m + n*p + 1, 0] && IntBinomialQ[a, b, c, n, m, p, x]

Rule 327

Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a + b*x^n
)^(p + 1)/(b*(m + n*p + 1))), x] - Dist[a*c^n*((m - n + 1)/(b*(m + n*p + 1))), Int[(c*x)^(m - n)*(a + b*x^n)^p
, x], x] /; FreeQ[{a, b, c, p}, x] && IGtQ[n, 0] && GtQ[m, n - 1] && NeQ[m + n*p + 1, 0] && IntBinomialQ[a, b,
 c, n, m, p, x]

Rubi steps \begin{align*} \text {integral}& = \frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {1}{4} a \int \frac {x^8}{\sqrt [4]{a+b x^4}} \, dx \\ & = \frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}-\frac {\left (5 a^2\right ) \int \frac {x^4}{\sqrt [4]{a+b x^4}} \, dx}{32 b} \\ & = -\frac {5 a^2 x \left (a+b x^4\right )^{3/4}}{128 b^2}+\frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {\left (5 a^3\right ) \int \frac {1}{\sqrt [4]{a+b x^4}} \, dx}{128 b^2} \\ & = -\frac {5 a^2 x \left (a+b x^4\right )^{3/4}}{128 b^2}+\frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {\left (5 a^3\right ) \text {Subst}\left (\int \frac {1}{1-b x^4} \, dx,x,\frac {x}{\sqrt [4]{a+b x^4}}\right )}{128 b^2} \\ & = -\frac {5 a^2 x \left (a+b x^4\right )^{3/4}}{128 b^2}+\frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {\left (5 a^3\right ) \text {Subst}\left (\int \frac {1}{1-\sqrt {b} x^2} \, dx,x,\frac {x}{\sqrt [4]{a+b x^4}}\right )}{256 b^2}+\frac {\left (5 a^3\right ) \text {Subst}\left (\int \frac {1}{1+\sqrt {b} x^2} \, dx,x,\frac {x}{\sqrt [4]{a+b x^4}}\right )}{256 b^2} \\ & = -\frac {5 a^2 x \left (a+b x^4\right )^{3/4}}{128 b^2}+\frac {a x^5 \left (a+b x^4\right )^{3/4}}{32 b}+\frac {1}{12} x^9 \left (a+b x^4\right )^{3/4}+\frac {5 a^3 \tan ^{-1}\left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{256 b^{9/4}}+\frac {5 a^3 \tanh ^{-1}\left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{256 b^{9/4}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.52 (sec) , antiderivative size = 98, normalized size of antiderivative = 0.78 \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=\frac {2 \sqrt [4]{b} x \left (a+b x^4\right )^{3/4} \left (-15 a^2+12 a b x^4+32 b^2 x^8\right )+15 a^3 \arctan \left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )+15 a^3 \text {arctanh}\left (\frac {\sqrt [4]{b} x}{\sqrt [4]{a+b x^4}}\right )}{768 b^{9/4}} \]

[In]

Integrate[x^8*(a + b*x^4)^(3/4),x]

[Out]

(2*b^(1/4)*x*(a + b*x^4)^(3/4)*(-15*a^2 + 12*a*b*x^4 + 32*b^2*x^8) + 15*a^3*ArcTan[(b^(1/4)*x)/(a + b*x^4)^(1/
4)] + 15*a^3*ArcTanh[(b^(1/4)*x)/(a + b*x^4)^(1/4)])/(768*b^(9/4))

Maple [A] (verified)

Time = 4.37 (sec) , antiderivative size = 126, normalized size of antiderivative = 1.01

method result size
pseudoelliptic \(\frac {128 b^{\frac {9}{4}} \left (b \,x^{4}+a \right )^{\frac {3}{4}} x^{9}+48 a \,b^{\frac {5}{4}} x^{5} \left (b \,x^{4}+a \right )^{\frac {3}{4}}-60 a^{2} x \,b^{\frac {1}{4}} \left (b \,x^{4}+a \right )^{\frac {3}{4}}-30 \arctan \left (\frac {\left (b \,x^{4}+a \right )^{\frac {1}{4}}}{b^{\frac {1}{4}} x}\right ) a^{3}+15 \ln \left (\frac {-b^{\frac {1}{4}} x -\left (b \,x^{4}+a \right )^{\frac {1}{4}}}{b^{\frac {1}{4}} x -\left (b \,x^{4}+a \right )^{\frac {1}{4}}}\right ) a^{3}}{1536 b^{\frac {9}{4}}}\) \(126\)

[In]

int(x^8*(b*x^4+a)^(3/4),x,method=_RETURNVERBOSE)

[Out]

1/1536*(128*b^(9/4)*(b*x^4+a)^(3/4)*x^9+48*a*b^(5/4)*x^5*(b*x^4+a)^(3/4)-60*a^2*x*b^(1/4)*(b*x^4+a)^(3/4)-30*a
rctan(1/b^(1/4)/x*(b*x^4+a)^(1/4))*a^3+15*ln((-b^(1/4)*x-(b*x^4+a)^(1/4))/(b^(1/4)*x-(b*x^4+a)^(1/4)))*a^3)/b^
(9/4)

Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.28 (sec) , antiderivative size = 234, normalized size of antiderivative = 1.87 \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=\frac {15 \, \left (\frac {a^{12}}{b^{9}}\right )^{\frac {1}{4}} b^{2} \log \left (\frac {125 \, {\left ({\left (b x^{4} + a\right )}^{\frac {1}{4}} a^{9} + \left (\frac {a^{12}}{b^{9}}\right )^{\frac {3}{4}} b^{7} x\right )}}{x}\right ) - 15 i \, \left (\frac {a^{12}}{b^{9}}\right )^{\frac {1}{4}} b^{2} \log \left (\frac {125 \, {\left ({\left (b x^{4} + a\right )}^{\frac {1}{4}} a^{9} + i \, \left (\frac {a^{12}}{b^{9}}\right )^{\frac {3}{4}} b^{7} x\right )}}{x}\right ) + 15 i \, \left (\frac {a^{12}}{b^{9}}\right )^{\frac {1}{4}} b^{2} \log \left (\frac {125 \, {\left ({\left (b x^{4} + a\right )}^{\frac {1}{4}} a^{9} - i \, \left (\frac {a^{12}}{b^{9}}\right )^{\frac {3}{4}} b^{7} x\right )}}{x}\right ) - 15 \, \left (\frac {a^{12}}{b^{9}}\right )^{\frac {1}{4}} b^{2} \log \left (\frac {125 \, {\left ({\left (b x^{4} + a\right )}^{\frac {1}{4}} a^{9} - \left (\frac {a^{12}}{b^{9}}\right )^{\frac {3}{4}} b^{7} x\right )}}{x}\right ) + 4 \, {\left (32 \, b^{2} x^{9} + 12 \, a b x^{5} - 15 \, a^{2} x\right )} {\left (b x^{4} + a\right )}^{\frac {3}{4}}}{1536 \, b^{2}} \]

[In]

integrate(x^8*(b*x^4+a)^(3/4),x, algorithm="fricas")

[Out]

1/1536*(15*(a^12/b^9)^(1/4)*b^2*log(125*((b*x^4 + a)^(1/4)*a^9 + (a^12/b^9)^(3/4)*b^7*x)/x) - 15*I*(a^12/b^9)^
(1/4)*b^2*log(125*((b*x^4 + a)^(1/4)*a^9 + I*(a^12/b^9)^(3/4)*b^7*x)/x) + 15*I*(a^12/b^9)^(1/4)*b^2*log(125*((
b*x^4 + a)^(1/4)*a^9 - I*(a^12/b^9)^(3/4)*b^7*x)/x) - 15*(a^12/b^9)^(1/4)*b^2*log(125*((b*x^4 + a)^(1/4)*a^9 -
 (a^12/b^9)^(3/4)*b^7*x)/x) + 4*(32*b^2*x^9 + 12*a*b*x^5 - 15*a^2*x)*(b*x^4 + a)^(3/4))/b^2

Sympy [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 3.62 (sec) , antiderivative size = 39, normalized size of antiderivative = 0.31 \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=\frac {a^{\frac {3}{4}} x^{9} \Gamma \left (\frac {9}{4}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {3}{4}, \frac {9}{4} \\ \frac {13}{4} \end {matrix}\middle | {\frac {b x^{4} e^{i \pi }}{a}} \right )}}{4 \Gamma \left (\frac {13}{4}\right )} \]

[In]

integrate(x**8*(b*x**4+a)**(3/4),x)

[Out]

a**(3/4)*x**9*gamma(9/4)*hyper((-3/4, 9/4), (13/4,), b*x**4*exp_polar(I*pi)/a)/(4*gamma(13/4))

Maxima [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 189, normalized size of antiderivative = 1.51 \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=-\frac {5 \, a^{3} {\left (\frac {2 \, \arctan \left (\frac {{\left (b x^{4} + a\right )}^{\frac {1}{4}}}{b^{\frac {1}{4}} x}\right )}{b^{\frac {1}{4}}} + \frac {\log \left (-\frac {b^{\frac {1}{4}} - \frac {{\left (b x^{4} + a\right )}^{\frac {1}{4}}}{x}}{b^{\frac {1}{4}} + \frac {{\left (b x^{4} + a\right )}^{\frac {1}{4}}}{x}}\right )}{b^{\frac {1}{4}}}\right )}}{512 \, b^{2}} - \frac {\frac {5 \, {\left (b x^{4} + a\right )}^{\frac {3}{4}} a^{3} b^{2}}{x^{3}} + \frac {42 \, {\left (b x^{4} + a\right )}^{\frac {7}{4}} a^{3} b}{x^{7}} - \frac {15 \, {\left (b x^{4} + a\right )}^{\frac {11}{4}} a^{3}}{x^{11}}}{384 \, {\left (b^{5} - \frac {3 \, {\left (b x^{4} + a\right )} b^{4}}{x^{4}} + \frac {3 \, {\left (b x^{4} + a\right )}^{2} b^{3}}{x^{8}} - \frac {{\left (b x^{4} + a\right )}^{3} b^{2}}{x^{12}}\right )}} \]

[In]

integrate(x^8*(b*x^4+a)^(3/4),x, algorithm="maxima")

[Out]

-5/512*a^3*(2*arctan((b*x^4 + a)^(1/4)/(b^(1/4)*x))/b^(1/4) + log(-(b^(1/4) - (b*x^4 + a)^(1/4)/x)/(b^(1/4) +
(b*x^4 + a)^(1/4)/x))/b^(1/4))/b^2 - 1/384*(5*(b*x^4 + a)^(3/4)*a^3*b^2/x^3 + 42*(b*x^4 + a)^(7/4)*a^3*b/x^7 -
 15*(b*x^4 + a)^(11/4)*a^3/x^11)/(b^5 - 3*(b*x^4 + a)*b^4/x^4 + 3*(b*x^4 + a)^2*b^3/x^8 - (b*x^4 + a)^3*b^2/x^
12)

Giac [F]

\[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=\int { {\left (b x^{4} + a\right )}^{\frac {3}{4}} x^{8} \,d x } \]

[In]

integrate(x^8*(b*x^4+a)^(3/4),x, algorithm="giac")

[Out]

integrate((b*x^4 + a)^(3/4)*x^8, x)

Mupad [F(-1)]

Timed out. \[ \int x^8 \left (a+b x^4\right )^{3/4} \, dx=\int x^8\,{\left (b\,x^4+a\right )}^{3/4} \,d x \]

[In]

int(x^8*(a + b*x^4)^(3/4),x)

[Out]

int(x^8*(a + b*x^4)^(3/4), x)